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Continuous measurement with traveling wave probes 
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We consider the use of a traveling wave probe to continuously measure the quantum state of an 
atom in free space. Unlike the more familiar cavity QED geometry, the traveling wave is intrinsi- 
cally a multimode problem. Using an appropriate modal decomposition we determine the effective 
measurement strength for different atom-field interactions and different initial states of the field. 
These include the interaction of a coherent-state pulse with an atom, the interaction of a Fock-state 
pulse with an atom, and the use of Faraday rotation of a polarized laser probe to perform a QND 
measurement on an atomic spin. 
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I. INTRODUCTION 

Current research to create quantum information pro- 
cessors has forced a reexamination of the underlying de- 
scription of these devices. In order for us to gain in- 
formation about quantum systems they must be mea- 
sured. Whereas the standard picture of such measure- 
ments involves the "collapse of the wave function" follow- 
ing Von Neumann's projection postulate, such strongly 
measuring probes are rarely implemented in the labora- 
tory. More typically a continuous probe interacts with 
the system which is then detected as a macroscopic sig- 
nal. Examples include the probing of a quantum dot with 
a single electron transistor (SET) and the measure- 
ment of the position of a micro- mechanical cantilever by 
monitoring the modulation of a reflected laser beam [2( • 
The formalism of quantum mechanics provides a num- 
ber of different approaches for analyzing such situations. 
Scattering theory employs Green's function input-output 
relations to describe the evolution of the probe asymp- 
totically, both before and after its interaction with the 
system under examination. Alternatively, the theory of 
quantum trajectories provides a dynamical description 
of the quantum system being measured, conditioned on 
the continuous information being collected via the probe. 
As laboratory developments give us access to the control 
and manipulation of quantum systems these descriptions 
become ever more relevant. The quantum trajectory ap- 
proach has the advantage of directly tying the dynamics 
of the system's evolution to the measurement record. The 
ability to do this is essential when implementing adaptive 
measurement and control strategies employing feedback 



The standard paradigm for continuous measurement is 
cavity QED. The dynamics of a cavity mode of the elec- 
tromagnetic field (perhaps coupled to an atom) are mon- 
itored by a partially transmitting mirror 0, Q ■ Input- 
output scattering theory, suited specifically to the lan- 
guage of optical elements ^J, is used to connect the 
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intracavity dynamics with those of the traveling signal. 
In order to translate the typically discrete information of 
the individually transmitted photons into continuous in- 
formation one considers a homodyne or heterodyne mea- 
surement, in which the signal is mixed with a macroscopic 
local oscillator. The result is a stochastic Schrddingcr 
equation which describes both the localization of the 
quantum state conditioned on the measurement and also 
the effect of quantum "back-action noise" 

In many applications, there is no confining cavity to 
set the interaction strength between the system under 
observation and the probe. Though a "quantization vol- 
ume" may be employed as a calculational tool, its imag- 
ined characteristics should not determine physically rel- 
evant quantities. For example, Milburn et al. 2] mod- 
elled continuous observation of a moving cantilever mon- 
itored by the modulation of a reflected laser beam. They 
considered the cantilever to be a mirror, which, part- 
nered with an imagined partially transmitting surface, 
formed a leaky optical cavity. Under the assumption that 
the transmission rate of light through the fictitious mir- 
ror is much faster than the characteristic rate at which 
the cantilever moves, the cavity could be adiabatically 
eliminated from the dynamics. This led to a stochas- 
tic Schrodinger equation for the continuously observed 
cantilever alone. The unphysical quantization volume, 
however, still appeared implicitly in this equation. 

The key parameter that characterizes the dynamics 
of a continuously observed system is the "measurement 
strength" , which we will denote n. It determines the 
rate at which information is gathered about the system 
and consequently sets the scale at which effects such as 
quantum back-action become significant. For example, 
Bhattacharya et al. ^| have shown that for sufficiently 
macroscopic systems there is a window of values for k 
such that continuous observation can localize the proba- 
bility distribution to a quantum trajectory that faithfully 
tracks the classically predicted trajectory, with minimal 
quantum noise. Another example is the continuous mea- 
surement of ensembles of atoms, controlled through their 
collective interaction with a common probe, to produce 
nonclassical spin squeezed states 0, hj. These effects 
depend crucially on k and its relation to the other rates 



2 



governing the system dynamics. 

In this article we formulate the problem of continu- 
ous measurement by a traveling wave probe. We derive 
a master equation describing the situation in which the 
system is monitored by the probe, but the measurement 
result is not recorded. This allows us to to identify the 
important characteristic scales of the problem without 
reference to a particular measurement scheme. We be- 
gin in Sec. m by establishing the necessary formalism for 
treating propagating fields in quantum optics, in con- 
trast to the more familiar closed cavity problems. We 
apply this formalism in Sec. lIIII to the fundamental prob- 
lem of a two-level atom interacting with a resonant laser 
field. When the field is treated classically this leads to 
Rabi flopping, but when treated quantum mechanically 
the laser not only manipulates the atom but also acts 
to continuously measure it. We determine the rate at 
which the measurement back-action leads to decoherence 
in the atomic system. This is contrasted with the evo- 
lution when the atom is coupled to a resonant electro- 
magentic pulse with a fixed photon number n. In partic- 
ular, we explore the circumstances under which we can 
recover the usual Jaynes-Cummings solution for a two- 
level atom coupled to a single mode and show how 
the behavior for a single photon diverges from this solu- 
tion. In Sec. HH we consider continuous measurement of 
an atomic spin through the Faraday rotation of an off- 
resonant laser field. This process has been used to gen- 
erate spin-squeezed states in atomic ensembles |0, Tl6| . 
We conclude and summarize our results in Sec. 



II. QUANTUM DESCRIPTION OF 
PROPAGATING FIELDS 

Classically, when considering quasimonochromatic 
propagating fields, it is natural to model the evolution 
of the system as a function of the propagation direction, 
z. The field at z can then be decomposed into a com- 
plete set of orthonormal temporal modes which act lo- 
cally. One might be tempted to describe the quantum 
fields in an analogous manner by quantizing the tem- 
poral modes at fixed position, [a(t), a*(tf)] = S(t — t'). 
The field operator could then be decomposed into a com- 
plete set of orthonormal mode functions, 4>i(t), so that 
a(t) = Y^i <f>i{t) Cj; with [ci,cj] = Sij. Boundary condi- 
tions at some initial plane could then used to restrict the 
mode content, possibly to a single temporal mode. 

This approach was taken by van Enk and Kimble |l7| 
and also by Gea-Banacloche 18] who considered an anal- 
ogous problem to the one we address here . They studied 
how errors were generated in quantum logic operations 
due to the fact that control pulses are not truly classical 
and can become entangled with the atoms with which 
they are interacting. Their analysis led to an effective 
single temporal mode theory. Though some of their con- 
clusions are correct, one must take great care to under- 
stand the regimes under which this formalism is applica- 



ble. Consider, for example, a single photon pulse inter- 
acting with a localized two-level atom. Let us suppose 
that the duration of the pulse is short compared to the 
natural lifetime of the atom in its excited state but suf- 
ficiently long to be considered quasimonochromatic and 
on resonance. Defining creation and annihilation oper- 
ators for the single temporal mode associated with this 
pulse, the Hamiltonian, in the rotating wave approxi- 
mation, appears to have the familiar Jaynes-Cummings 
form, 

H = hg (a<r + + a^(j_) , (1) 

where o± are the usual raising and lowering operators 
associated with the two-state atom. Given the atom ini- 
tially in its ground state, the solution leads to quantum 
Rabi oscillations, 

\m) = cos( 5 t) |s) |1) -isin(gt) \e) |0) . (2) 

This falsely predicts the possibility of a single photon 
27r-pulse in free space, whereby the photon is perfectly 
absorbed and then reemitted into the original mode. In 
reality once the atom has absorbed the photon it will 
reemit into a mode consistent with its radiation pattern, 
not into the initial packet mode. That is, the single 
photon will be scattered. This emission must also obey 
causality; no information about the emitted photon can 
register on a distant detector at a space-like separated 
point. In the solution above, however, the atom both ab- 
sorbs from and emits into a spatially delocalized photon 
mode in free space, violating causality. 

The problems with causality arise from the faulty 
quantization procedure outlined above. Quantum fields 
must be defined over all space at an initial time (more 
generally on an initial space- like surface) . Unitarity then 
ensures that equal-time, not equal-space, commutation re- 
lations are preserved. Nonequal-time commutation re- 
lations cannot generally satisfy the canonical commuta- 
tion relations, being inconsistent with Poincare invari- 
ance [l^. The exception is for free fields, or fields that 
behave like them (e.g. fields traveling through matter 
whose response is approximately linear, or asymptotic 
"in" and "out" fields as used in scattering theory). 

We review here a formalism appropriate for tre ating 
the quantum optics of paraxial propagating fields |20| . 
Consider a quasimonochromatic paraxial beam with fre- 
quency loq and wave number ko. We write the positive 
frequency component, 

£*+>(M) = s(z,t)M^yy ikoz ~ Uot) , (3) 

where exp[i(koz — uJot)] is the "carrier wave", (j)T(x,y) 
is the "transverse mode" (e.g. Gaussian), and £(z,t) 
is the slowly varying envelope, meaning its spatio- 
temporal variation is much slower than the carrier wave- 
length/frequency. We have ignored both diffraction and 
the vector nature of the field. It is easy to show that the 
free space wave equation becomes, 

(l+4) £ <^> = ° (4) 
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for the envelope, whose solution is E (z,t) = £ (z — ct, 0), 
i.e. propagation of the pulse envelope. We quantize by 
replacing the field envelope with a scaled operator, 



£(z,t) => v/27r^ *(z,i), 



(5) 



which satisfies the canonical equal-time commutation re- 
lation m, 



V(z,t),W{z',t) =5{z-z') 



(6) 



This commutator is equivalent to that of a nonrelativistic 
massive Bose gas in one dimension. Here the carrier wave 
plays the role of the rest mass and the slowly varying 
envelope plays the role of small fluctuations around the 
mass shell. The free field Hamiltonian (removing the 
carrier wave energy) takes the form, 



H F = c J dz¥(z) ^-ih 



(7) 



whose Heisenberg equation of motion gives the wave 
equation above. This Hamiltonian is nothing but the sec- 
ond quantized version of the energy of a photon E = cp. 

Consider an atom interacting with the field. In the 
electric-dipole and rotating wave approximation the in- 
teraction Hamiltonian is, 

Haf = J d 3 x |$(x)| 2 d(i;W(x) ( 7 + +M-\x)a- 

(8) 

where d is the dipole matrix element and $(x) is the 
atom's center of mass wave function. We take the atom 
to be trapped, having center of mass wave function 
<f>(x) = fr(x, y) j25|. Then the interaction Hamil- 

tonian becomes, 



H 



AF 



where J dxdy\f T (x,y)\ 2 (f> T {x,y) = -^=, (9) 



A being the effective area of the mode interacting with 
the atom. Let us go to the interaction picture by in- 
cluding the free evolution of the atom and field in the 
interaction Hamiltonian. Assuming the carrier wave is 
on resonance, 



H AF (t) = d 



A 



dz 



\.f L {z)\ 2 (y{z-ct)a + (10) 



+$ + (z - ct)a- 



Finally, given a set of orthonormal functions (" longitu- 
dinal modes") {<j>i(z)}, chosen to be real without loss of 
generality, 



(11) 



(f>i(z - ct) ( 



a,ia + + a- a- 



(12) 



The Hamiltonian in Eq. (|12[) describes the interaction 
of each longitudinal mode as it propagates past the 
atom(Fig.^|. Assuming this time scale is much shorter 
than any other dynamical scale in the problem, it is ap- 
propriate to make the Markov approximation, whereby 
one coarse grains over the time scale over which the sys- 
tem and reservoir retain "memory" . To this end, we 
break up the z-axis into slices of size Az, each the ex- 
tent of the atom wave packet (e.g. the rms of the prob- 
ability density). The Markov approximation will hold 
if the transit time of the field across the atomic packet, 
At = Az/c, is much smaller than any time scale over 
which the atom changes. This is certainly an excellent 
approximation. Since we will not consider dynamics on 
a time scale smaller than At, we can approximate the set 
of atomic wavepackets, centered at each coarse grained 
slice, as a complete orthonormal set. That is, each slice 
is an approximate delta function. Then normalization of 
both the atomic wave packet and mode functions com- 
bine to give 



dz \/l(z)\ 4>i{z - ct) 



1 



=9i(t) 



Az 

(i - l)At < t < iAt 
otherwise 



(13a) 
(13b) 



Under this approximation the Hamiltonian takes the 
form, 



HAF{t) = Hg^2@i(t) (a.i& + + ajo-_ 



(14) 



where hg = dE v 



2ttHujo 
AcAt 



This result has a clear interpretation. The traveling wave 
configuration is effectively multimode, with each mem- 
ber of the set being a traveling packet "mode-matched" 
to the atom. The coupling constant g depends on this 
mode volume. This picture is equivalent to a model of 
decoherence discussed by Brun 26] in which a "flying 
qubit" passes over a "system qubit" , the former acting 
as an irreversible reservoir (through its continuous spatial 
degrees of freedom) to carry information away from the 
system, thereby leading to decoherence. In our problem 
a given harmonic oscillator (mode of the electromagnetic 
field) flies over the qubit, becomes entangled with it, and 
then flies away. This too leads to decoherence, as we 
describe in the next section. 
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FIG. 1: An initial traveling wave pulse (solid line) is broken 
up into many smaller modes {(f>i(z)} (dashed line). The wave 
packet of the interacting atom (gray) has the same width as 
the mode functions (Az). 



III. THE RABI INTERACTION FOR 
TRAVELING WAVES 

A. Interaction with a laser beam 

We consider first the case of a resonant laser beam 
interacting with a trapped two-level atom. The state 
of the field is described by a tensor product of identical 
coherent states for each traveling mode packet, 



/ beam 



(15) 



where the amplitude is given by the mean number of 
photons in time slice At, \a\ 2 = PAt/(hujo), with P the 
power of the laser beam. More precisely, the state of the 
beam is a statistical mixture of states of the form above, 
averaged over the common, but unknown, phase of the 
complex amplitude a, as described by van Enk and Fuchs 
|2l| . The actual value of the phase plays no role in the 
analysis to follow, so we choose it to be fixed with no loss 
of generality. In order to distinguish coherent effects from 
decoherence we transform by a unitary displacement of 
the field states to the vacuum, 

^({a}) |¥> , A^D-\{a})AD({a}). (16) 
In this picture, 

HAF(t) = hga (<7+ + <x_) + hg^2 <di(t) [a l <j + + a\a^j 

i 

= -Hcoh + HaV- (17) 

The coherent term is classical Rabi flopping at a fre- 
quency Q, = 2ga = dW8%I/Wc, with / the beam inten- 
sity (cgs units, as used throughout). The second term is 



the atom- vacuum coupling for the travelin g^ wa ve modes 
only (i.e. the paraxial modes of the beam) |22fl . 

We can now proceed with the standard Markov anal- 
ysis to derive the Master equation. The initial atom- 
vacuum state is uncorrelated. After a time At, one of 
the modes becomes entangled with the atom through the 
atom- vacuum coupling. The Linblad ("jump") operator, 
L, is defined by |2q . 



(1| U AV (At) |0) = —At (1| H AV |0) = Lv At. (18) 
Plugging in the Hamiltonian from Eq. p7jl. we arrive at, 



2-Khulg - — 
A „ U _ 



(19) 



where k = d 2 (^l) = T 



r 



. A I 



is the measurement strength, T the spontaneous emis- 
sion rate and cr e ff the effective cross section for scattering 
out of the paraxial modes. By determining these jump 
operators, we compactly derive the master equation, 
equivalent to that obtained using a the usual system- 
reservoir approach after tracing over the unmeasured 
bath dHHH, 

= -r [#coh,/3] - | {d+o-_, p] + na_pa + . (20) 

This equation has a familiar and appealing form. It is 
none other than the master equation for a decaying laser- 
driven atom 0|, but with r — > k. The decay rate is 
due to the entanglement between the atom and the laser 
modes. Note, k is independent of At, which acts as a 
fictitious "quantization volume" and so must be absent 
from any physical quantities such as the measurement 
strength. 

The measurement strength is also independent of the 
laser power \a\ 2 . In particular we may turn off the laser 
(a — > 0), and the measurement strength will remain the 
same. The ratio k/T — <t c h/A may thus be interpreted 
as the fraction of spontaneous emission into the parax- 
ial modes. In support of this interpretation note that 
the mode area A can never be made smaller than the 
diffraction limit A ~ 1/feg, so at most n ~ T. Moreover, 
once the beam becomes focused to such a small spot size, 
one can no longer neglect the vector nature of the atom 
field coupling which further decreases the measurement 
strength |27j . 

From Eq. I|19f) we can determine how continuous mea- 
surement by the laser beam acts to decohere the atom. 
For a paraxial beam we require that a e s/A 3> 1, so that 
diffraction effects are minor. Then decay due to entan- 
glement with the laser modes is small compared to decay 
due to spontaneous emission into 4-7T steradians. In agree- 
ment with the conclusions of 17. 18], errors in coherent 
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control pulses due to the quantum nature of the inter- 
action can be neglected if spontaneous emission is also 
negligible during the duration of the interaction. 



B. Interaction with a single photon 

In Sec. m we showed how a quantization procedure in 
terms of nonequal-time commutators can lead to a false 
prediction of single photon coherent Rabi flopping in free 
space. In this subsection we use our formalism to show 
how a quasimonochromatic and paraxial propagating sin- 
gle photon wave packet drives a two-level atom. 

We take the initial state of the system to be a single 
photon wave packet with the atom in its ground state, 



9{0))=V[f] |vac)®|<?) A . 



(21) 



The operator a'[f] = J dz^ (z) f (z) creates a delocalized 
single photon state with slowly varying pulse envelope 
f(z) [29j . For simplicity we use a square pulse of duration 
r = NAt. In this case we can expand f{z) in a symmetric 
sum of coarse-grained modes each having length Az = 
At/c, 



tm - 



an/] 



1 



jV-1 



(22) 



The state will evolve according to the Hamiltonian in 
Eq. I|14l) which commutes with the total number of ex- 
citations in the system. Neglecting, for now, the possi- 
bility of spontaneous emission into other transverse field 
modes, the total number of excitations will be preserved. 
The state at all times must then have the form 



(23) 



times t k < t j: shows that A k {t k -i) = A k (0) = N 1/2 . 
Inserting this result into Eq. H25cfl . 



A e (t k ) = A e {tk-x)c- -=. 

v N 



This equation admits a simple series solution, 



A e (t k ) = - 



is 1 — cr 



N 1-c 



(26) 



(27) 



Assuming viV > 1, i.e. the envelope f(z) is much 
broader than the coarse graining, we may take the limit 
of A e (t) as N -> oo holding r = NAt. This yields, 



A e (t) 



2i 



1 - e 



-Kt/2 



(28) 



The solution given in Eq. I|28|) is based on the funda- 
mental assumption that the evolution of the state is uni- 
tary, i.e. we consider a closed quantum system consisting 
of the atom and paraxial field modes. In the continuum 
limit, this yields an effective decay due to emission into 
the included paraxial modes at rate k, but it excludes de- 
cay into all others modes which, taken together, give a to- 
tal spontaneous emission rate T. Since we showed in Sec. 
(|III All that k <C r, this solution is not self-consistent. To 
rectify this, during the time interval At we must allow 
for a small probability of spontaneous emission, P S pont, 
into non-paraxial modes. By not including these modes 
in our system, the state ket evolves according to an 
effective non-Hcrmitian Hamiltonian [12J with decaying 
norm, = 1 — P spon . Eq. I|25c|) then reads, 



A e (t k ) = A e (t k - X )c& 



-7Ai/2 



iA k (t k -l)s, 



(29) 



where 7 is the spontaneous emission into all non-paraxial 
modes. Employing the initial condition and taking the 
limit At -> 0, NAt = r yields, 



Consider the evolution of the system over the short 
interval (t k ,t k + At], where t k = kAt,0 < k < N - 1. 
We can define a map for the state between two successive 
time steps, 



_ -iH k At/h 



M*k-l)> 



(24) 



where H k — hg (^a k a + + a\,<T- \ . Using ansatz (|23f) we 
are led to the recursion relations, 



A 3 ^ k (t k ) = Aj(t k -i), 

A k (t k ) = A k {t k ^i)c - iA e (t k -i)s, 
A e (t k ) = A e (t k -i)c - iA k (t k -i)s. 



(25a) 
(25b) 
(25c) 



Here s = sin V nAt and c = cos V nAt. The measurement 
strength k — g\/~Ai is the same as in Eq. (|19|l . 

These coupled algebraic equations can be solved for the 
amplitudes. Repeated application of the Eq. I|25a|) at all 



dt A ^ 



•-(7 + /t)A»(t) 



(30) 



Solving this equation with 7 = 0, using the initial condi- 
tion ^4 e (0) = 0, will give the same result as in Eq. I|28l) . 

Before considering the solution to this differential 
equation, consider the slightly altered situation in which 
the field starts in the vacuum state, and the atom in the 
excited state. In this case it is easy to see that the second 
term on the right side of Eq. H30JI will vanish. Then the 
equation becomes, 



dt 



A e (t) 



1 



(j + K)A e (t), 



(31) 



with the initial condition A e (0) = 1. This must give the 
standard exponential decay due to spontaneous emission 
in the vacuum 



A e {t) 



(32) 
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FIG. 2: The probability for a single photon in free space to 
excite a two-level atom (solid line) differs significantly from 
the standard Rabi flopping solution, in which an atom is cou- 
pled to a photon in a single-mode cavity ( dashed line). The 
parameters used here are Tt = 2.5, and k/T — 1/50 (see 
text). 

This allows us to equate r = 7 + k, where k is again seen 
to be the spontaneous emission rate due to the contribu- 
tion of the modes in the paraxial beam. 
The general solution to Eq. i|3TJ|) is, 

A e (t)=-^(l-e-™/*), (33) 

holding for times < t < r. The probability of the atom 
being in the excited state is then, 

P e (t) = 4^ (l - e- rt/2 ) 2 , (34) 

during the same interval. This is a monotonically increas- 
ing function of t, and so achieves its maximum at the 
upper limit t — r, after which the excitation probability 
can only decay (Fig. Then the maximum probability 
for any r can be found by solving, 

^-Pe(r) = 0->Tt^ 2.5. (35) 

At this point the probability is P e ~ .8re/r, which is 
necessarily less than 1 given that k/T < 1 as previously 
discussed. Thus, in the paraxial approxiamtion, no single 
photon pulse can be constructed that excites an atom 
with certainty. 

This result is, of course, not surprising. Symmetry en- 
sures that the only single photon pulse capable of exciting 
an atom with unit probability is the time reversal of a 
spontaneously emitted packet [3(| • Such a pulse repre- 
sents the vector spherical harmonic associated with the 
atomic radiation pattern (assuming emission on a tran- 
sition with well defined angular momentum), a field not 



captured in the scalar paraxial approximation. More- 
over, even including the vector nature of the field be- 
yond paraxial will not be sufficient to yield high excita- 
tion probability. The field must be well "mode matched" 
to the the atom's radiation pattern to give strong cou- 
pling between the atom and a single photon in free space 



C. Interaction with a large ?i-Photon Fock Pulse 

In the last subsection we showed that in free space the 
coupling of a single photon to an atom will not lead to 
coherent Rabi oscillations. In contrast, for an n-photon 
Fock state with a large n, we expect the system to be 
dominated by stimulated emission. A rigorous treat- 
ment of this probem using the Bethe-Ansatz was given 
by V. I. Rupasov and V. I. Yudson |28j. We show here 
how this phenomenon is recovered simply in the present 
formalism. 

In general, an arbitrary division may be envisioned in 
Hilbert space that separates system from reservoir. In 
Sec. [n] the system was chosen to be the paraxial field 
modes plus the atom. Alternatively we may take the 
system to consist of the atom interacting with the single 
pulse mode defined by creation operator fit [/] in Eq. 122fl . 
with all other modes making up the environment. Such 
a choice can always be made, and cannot change the 
physics. We may then choose to ignore all environmental 
modes if Ft <C 1, with the caveat that any effects that 
occur in the system must be have scales much larger than 
Tt in order to be considered valid. 

Consider a pulse of length r such that rr C 1. Fur- 
ther, take the state of the system to only have excitations 
in this pulse mode such that, 

|^(O)) = -L(at[/]) n 0|O, 5 ) = K ff ). (36) 

The rest of the modes are in the vacuum state and are 
treated as an environment. Ignoring terms of order Tt 
the system evolves under the single mode Hamiltonian 

& = fifleff(o[/]* + + fit [/]*_), (37) 

This is none other that the Jaynes-Cummings Hamil- 
tonian restricted to an initial manifold with n exci- 
tations and having effective coupling constant hg e s = 
d^j '2-KhjjQ /Vpuisc , where Vpulse = Act is the pulse vol- 
ume. The system undergoes the familiar coherent Rabi 
flopping within the two dimensional manifold, as in Eq. 
©, 

\ip{t)) = cos(g e sy/nt) \g) \n) - i sm.(g e sy/nt) |e) \n— 1} . 

(38) 

This solution applies to the single photon case as well. 
The probability amplitude from Eq. (|33|) limits to, 

A e (t)^-i9est(l + 0(Tt)), (39) 
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which is the correct limit of sinusoidal Rabi oscillation. 
However during the pulse duration not even a single os- 
cillation can occur since g e g = \/k/t -C 1/t . Thus, for 
true oscillations to occur, one must have 



1 .A 
lesyn — or nTr <. . 

~ T ~ fcff 



(40) 



The last inequality can be interpreted as saying that the 
mean number of photons emitted via stimulated emission 
into the pulse must dominate over spontaneous emis- 
sion, even when the spontaneous photons are paraxial. 
When these conditions hold we may consistently ignore 
all initially unoccupied modes, yet still recover dynamics 
in agreement with the usual Jaynes-Cummings Hamilto- 
nian. The multimode description of the field becomes 
superfluous. 

Finally, what is the measurement strength associated 
with probing an atom using a large photon number Fock 
state pulse? Unlike the coherent state case, the field does 
not factorize into uncorrelated temporal slices. In fact, 
when viewed in terms of the coarse grained modes, the 
Fock pulse in highly entangled. Detection of photons at 
the leading edge of the pulse will introduce new fluctua- 
tions in the trailing edge, which has yet to interact with 
the atom. This implies that the usual Markov approxi- 
mations do not hold. The Fock pulse is most naturally 
treated as part of the "system" rather than an "environ- 
ment" which continuously carries information away from 
the atom. 



IV. QND MEASUREMENT VIA FARADAY 
POLARIZATION SPECTROSCOPY 



The resonant interaction considered up to this point, 
though fundamental in nature, has little practical appli- 
cation to the problem of continuous measurement in free 
space since the measurement strength is always bounded 
from above by the spontaneous emission rate. We thus 
turn our attention to an off-resonant interaction. In par- 
ticular, we consider the problem of measuring a spin com- 
ponent of an atom through the Faraday effect wherein 
the linear polarization of a probe laser rotates by an 
amount proportional to the magnetization of the sample. 
This interaction provides for a "quantum nondemolition 
measurement" (QND) of the atom, i.e. the probe mea- 
sures an observable (here the atomic spin along the laser 
propagation direction) without perturbing that observ- 
able [3TI l32| . Such an interaction has been applied to 
ensembles of atoms to produce spin squeezed states 0] 
and to demonstrate entanglement between two spatially 
separated ensembles [l6|. 

We review here the basic physics associated with this 
measurement scheme. The interaction energy is that of 
an electric-polarizable particle in a vector field whose 
Hamiltonian may be written as, 



where a is the atomic polarizability tensor and E is the 
complex electric field amplitude. Expressing this equa- 
tion in terms of irreducible tensor components, the inter- 
action can be decomposed into an effective scalar, vector, 
and symmetric rank-2 contribution. We consider here 
the case of alkali atoms probed on the so-called D2 line 
S1/2 — ► -P3/2; f° r which the ground state atomic polariz- 
ability operator is [34| 



am 



(42) 



Here e± are the right and left helicity vectors relative to 
the quantization axis along the probe propagation direc- 
tion, a z = ||) (T I — I J.) (J-l is the Pauli spin operator for 
the ground state electron, and a\[ n is the atomic polariz- 
ability for fields with linear polarization. This equation 
has a simple interpretation. The first "scalar term" gives 
rise to an effect that depends solely on the field inten- 
sity whereas the term proportional to a z depends on the 
field elipticity. Clearl y th e irreducible rank-2 component 
vanishes here [33l l34i l35j | . 

Under the paraxial approximation we take the field 
to be approximately a plane wave with two polarization 
vectors. The quantum field associated with the complex 
amplitude is 



E (+) = J-^(a_e_ 



(43) 



where V is the effective quantization volume for the prop- 
agating mode. Substituting this into Eq. (|41|l . the auan- 
tum Hamiltonian becomes, 



Hi n t = 



27rai i: 
V 



■huj 



(44) 

where N± is the number operator for photons in the pos- 
itive or negative helicity states. The scalar term gives 
rise to an overall phase shift (index of refraction) and 
thus can be absorbed into the free field Hamiltonian. The 
vector term gives rise to the Faraday effect. Recognizing 
J z = (N + — N_)/2 as the total field helicity, the effective 
interaction Hamiltonian takes the QND form, 



hit 



27roui 
V 



(45) 



Under this Hamiltonian the photon spin becomes cor- 
related with the atom's magnetic moment and thus the 
laser polarization may act as a meter for the atomic spin. 
Since a z also commutes with the system-meter Hamilto- 
nian it is clearly a QND variable jlj|. 

Our Hamiltonian still contains the undefined quanti- 
zation volume V . To rectify this we follow the formalism 
introduced in Sec. [H] Introducing propagating modes 
of duration At so that V — > AcAt the Hamiltonian be- 



comes, 



(41) 



£ 2 ge 4 (i)i 2 a z 



(46) 
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where x = not\i n u/{cA) w (<7 /A)[r/(-2A)]. Here 
A is the (far)detuning from the atomic resonance with 
linewidth T, and ctq is the on resonance absorption cross- 
section for linear polarization. 

With the Hamiltonian so defined, the evolution of the 
system may be calculated. The system shall consist of 
an atom interacting with a laser beam that is linearly 
polarized in the x-direction. The initial state of the laser 
is then, 



l $ )probc = (l«> to 1% 



(47) 



where x and y label the two orthogonal linear polariza- 
tion modes. As in Sec. IjHI A|) |a| 2 = PAt/hu , with P 
the power in the beam. Given this initial state for the 
field, we may define the marginal density operator for the 
atom alone by tracing over the field at time t. Since the 
Hamiltonian (|41|) only couples the atom to the fc'th field 
mode in time interval (t& , tfc + At] , the reduced atomic 
state evolves during this interval as, 



p(t k + At) = Tr fe Ui \a) ix (a\ lx ® |0> 



® p{th)U} 



(48) 

In the interaction picture the unitary evolution becomes, 



U t \a) lx \0) ly =e-^^-^\a) tx \0) ty , 

= \acos(xa z )) xi \-asin( X a z )) yl . (49) 

This expression should be interpreted in the context of 
a matrix element, given that atomic operators appear in 
the labels for the field kets. Clearly, the field and atomic 
states will generally become entangled by the interaction. 

The map in Eq. (|48|l constitutes a continuous measure- 
ment on the field. To see this explicitly, we must expand 
in powers of At, and since \a\ 2 oc At, this is equivalent 
to an expansion of the state kets in Eq. Ij49(l in powers of 
a. In the Fock (photon number) basis we have, 



Ui \a) ix 10), 



2 



JO)* |0) yi 

+ acos( X a,)|l)„|0), jl 

-asin(x<Tz) |0) xi \l) yl 

Substituting this back into Eq. I|48fl one finds, 

p{t k + At) = p(t k ) + \a\ 2 [-p(t k ) 

+ sin{xo z )p{tk) sin(x<5- z ) 

+ cos(x<7 z )p(t fc )cos(x<5- z )] 



(50) 



(51) 



This form may be further simplified since \ ^ 1 f° r a 
single atom probed by a far off resonance laser. Taking 
the limit, At — > 0, and keeping terms to second order in 
X gives the master equation, 



dp 
~dt 



= -2 [<Jz 



where the measurement strength can be easily identified 
from the familiar Linblad form of the mater equation |T^ 



hujQ ' 



(53) 



(52) 



Note that the steady state solutions to Eq. (|52J) are the 
Dicke states, so this technique provides a QND measure- 
ment of spin. This expression was also derived by Thom- 
sen and Wiseman in the context of control of atom-laser 
coherence 36]. 

The physics of the continuous measurement under con- 
sideration here differs substantially from that of the res- 
onant Jaynes-Cummings interaction studied in Sec. II I II 
Here the measurement strength depends upon the power 
in the beam, whereas the previous result had no such de- 
pendence. For the resonant could be explained 
as arising from those spontaneous photons emitted into 
paraxial modes. In contrast, the dependence of the QND 
measurement on the laser power indicates that the mea- 
surement strength is due to the coherent redistribution 
of photons between the polarization modes in a manner 
depending on the atomic state. 



V. SUMMARY AND DISCUSSION 

Continuous quantum measurement, once studied only 
in gedanken experime nts, can now be realized in labo- 
ratory applications || E| . Such applications utilize the 
ability to continuously gather information from a probe 
coupled to a quantum system. The system state then 
evolves under a stochastic master equation, character- 
ized by the measurement strength. In this article we 
examined how one may derive the measurement strength 
for paraxial laser probes in free space by considering two 
examples. 

We considered first the fundamental system of an atom 
coupled to a laser beam. Classically the atom-laser in- 
teraction leads to Rabi flopping which is often used to 
manipulate coherent superpositions of atomic states. In 
particular laser control pulses can be used to implement 
quantum logic gates [33. Quantum mechanically con- 
tinuous measurement of the atom by the quantum laser 
pulse can lead to entanglement between the system and 
the probe, inducing errors in logic gate operation. We 
found that this is not a concern as the measurement 
strength is bounded from above by the total sponta- 
neous emission rate and can thus be neglected for short 
interaction times. In fact, the measurement rate may 
be attributed to spontaneous emission into the paraxial 
modes, which will always be small compared to emission 
into 4-7T steradians with the appropriate vector field. We 
also studied the distinction between continuous observa- 
tion by a coherent beam and by a pulse with well defined 
photon number. In the latter case correlations between 
the quantum fluctuations at different points in the pulse 
disallow a clean separation between system and probe. 
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Instead, the entire atom-field system will coherently Rabi 
flop if the number of photons is sufficiently large. For a 
single photon, however, we showed that Rabi flopping 
cannot occur in free space. Such a result is essential lest 
a distant detector be able to instantly sense the pres- 
ence of the atom. This highlights the need to treat the 
propagation of quantum fluctuations in a traveling wave 
geometry with great care to avoid contradictions with 
causality. 

Another mechanism allowing for continuous measure- 
ment of an atom is the Faraday rotation of an off-resonant 
polarized beam, induced by the atom's magnetic mo- 
ment. This interaction produces a QND measurement 
of the atomic spin. The measurement is due to stimu- 
lated redistribution of photons between pairs of occupied 
modes, in contrast to the resonant Rabi interaction where 
measurement results from spontaneous emission. Since 
the measurement strength is proportional to the laser 
intensity it can be made much larger than the sponta- 
neous emission rate. Of course spontaneous emission is 
not the only source of noise. In order to see the effects 
of quantum back-action the quantum "projection" noise 
|38| must be large compared all other noise sources, such 
as those due to photo-detection. This implies that an ex- 
periment in which the back-action is important must en- 
hance the coupling constant without increasing the shot- 



noise, possibly through the use of an optical cavity or by 
maintaining a large ensemble of atoms collectively cou- 
pled to the probe. The latter approach has led to the 
observation of spin squeezing |l5j and ensemble entan- 
glement |l6j . In such a situation it would be interesting 
to understand the effective measurement strength by the 
probe on any one member of the ensemble. Exploring 
this question will require us to consider how quantum 
correlations are shared within a multipartite system in- 
teracting with a probe. We plan to address this question 
in future work. 
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